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Abstract 

It is shown that in the approximation of |p| = l the CP violation in B'^^ — > 
l^l^ decays vanishes in SM. In a 2HDM with CP violating phases and MSSM the 
CP asymmetries depend on the parameters of models and can be as large as 40% 
for B^ and 3% for B^. An observation of CP asymmetry in the decays would 
■ unambiguously signal the existence of new physics. 
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The flavor changing neutral current process, Bd,s — ^ (1=/^,''"), has attracted a lot 
^ : attention since it is very sensitive to the structure of SM and potential new physics beyond 
0\ , SM and was shown to be powerful to shed light on the existence of new physics before 
possible new particles are produced at colhders[|l], 0, 0]. In a very large region of parameter 
space supersymmetric(SUSY) contributions were shown to be easy to overwhelm the SM 
^ contribution ^, ^ and even reach, e.g., for l=/i, the experimental upper bound0 

5 ; Br{Bd 12+ fi-) < 6.8 X 10"^ {CL = 90%) 

Dh! Br{Bs < 2.0 X 10"^ {CL = 90%). (1) 

^ . In other words measuring the branching ratio of Bd^s 11 can give stringent con- 
straints on the parameter space of new models beyond SM, especially for that of the 
minimal supersymmetric standard model (MSSM) because of the tan^/3 dependence of 
^ \ SUSY contributions in some large tan/3 regions of the parameter space p|, |^. Comparing 
' with hadronic decays of B mesons, this process is very clean and the only nonperturbative 
quantity involved is the decay constant that can be calculated by using lattice, QCD sum 
rules etc. 

The results on CP violation in B^ - B^ mixing have been reported by the BaBar and 
Belle Collaborations in the ICHEP2000 Conference, which are consistent with the world 
average . More experiments on B physics have been planned in the present and future 
B factories P]. In the letter we study CP violation in Bd^s ^^^^ {^=f^y'^)y which might 
be measured in the near future. 

Obviously for the process Bd^s ^^^^ there are no direct CP violations since there are 
no strong phases in the decay amplitude 0, ^. But it is well known that CP violating 
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effects can survive after taking into account the mixing of tlie neutral mesons, and B^, 
in the absence of the strong phases. We will give a model-independent description for the 
CP violating effects of the process induced by mixing of 5° and and analyze them in 
SM and new models, a two Higgs doublet model(2HDM) with CP violating phases and 
MSSM. 

We need to know what kind of CP violating observables can be defined in the process. 
At first, direct CP violation, as noted above, is absent in this process. T-odd projection 
of polarization is a kind of useful tool to probe the CP violating effects, for example, in 
B ^ Xl^ l~ p|, |10|, [2^]. However for the process we are discussing here, we have actually 
only one independent momentum and one independent spin which can be chosen as those 
of Z", so no T-odd projections can be defined. Unlike the case generally discussed for 
hadronic final states, for example, that in Ref. |]TT|, the detected final states of and l~ 



of this process in experiments are basically two asymptotic energy-momentum eigenstates 
which are not CP eigenstates. Considering for instance B^ decays to in the rest 
frame of 5'', due to the energy-momentum conservation we denote the four-momenta of 
/~ and as p = {E,p) and p = {E,—p). Then the angular momentum conservation 
tells us that and final states are forbidden. Hence we are left with a pair of CP 
conjugated final states, and and two couple of corresponding CP conjugated 
processes. Therefore, we may define the time dependent CP asymmetries as 

Two corresponding time integrated CP asymmetries are 

^ /o" dt T{Bl^yAt) - U) - /o°° dt T{Bl,,yM - h) . ^ . , (A\ 
C dt T{Bl,^Xt) - h) + r dt nBl^^Xt) - /.) ' ' ^ ^ 

Where /i,2 = ^L,pL,Ry f is the CP conjugated state of /. 

The time evolutions of the initial pure B^ and B^ states are given by|l4 



\B',,yM)=9At)\B') + ^9-it)m, 



P 



\B',,yM) = -^9-it)\B') + g^it)\B'). (5) 



with g±{t) given by 



1 Am 
g+{t) = exp{--Tt - imt)cos{^t), 

1 Atti 
g~{t) = exp{—-Tt — imt)isin{—^t) (6) 

The absence of strong phases implies 

\Af\ = \Aj\, \Aj\ = \Af\ (7) 
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where Af{Af) =< f\neff\B^{B°) >. And the CPT invariance leads to 



Af _ (Af 



For simphcity, define 



Af \Aj 



Af _ ^j<t>f 1 ^ 



pe'^f^ ^ = xe'^-. (9) 
Af p 

From Eqs.(|]), (H), (||), (|^), and (||), it is straightforward to derive 

_ ^ AjB^t) ^ /) , _ |1 + x-'ptani^t)exp[t{-<Pf - 0. + f )]| 



Therefore, if 



or 



x^l (11) 



+ 0^ 7^ mod 2n7r, (12) 



( or equivalently /m(^^) 7^ 0, ) then r 7^ 1, i.e., one has CP violation. 

The effective Hamiltonian governing the process Bd^s —>■ has been given in Refs. 
0, 1^ . Using the effective Hamiltonian, we obtain by a straightforward calculation Q 

Af, ^ A, gqiy/l - 4mf + {Cq2 + 2m,Cio) 

where Aj = VthV^^ or Vj^VJ*, = mi/rriBO and Cj's are understood as Wilson coefficients 
at mB scale ^2], |1^, ^ |^. Because Cq.'s are proportional to m/ and Cio is independent 



of mi it follows from eq. ([I3|) that CP asymmetry in i^^.s /+Z is independent of the 
mass of the lepton. That is, it is the same for 1 = electron, muon and tau. 
In SM, one hasfig] 



P \Mu\ Xt 



(14) 



up to the correction smaller than or equal to order of 10^^, Cio is real, Cq^ = 0, and Cq, 
is negligibly small. So it follows from Eqs.([T3|), ( p!4D that a; = 1 and (j)f + (j)x=0. Therefore, 
there is no CP violation in SM. ^ If one includes the correction smaller than order of 10~^ 

^We have neglected the contributions, which is smaller than or equal to 10^"^ of the leading term, 
from the penguin diagrams with c and u quarks in the loop. It is true for both and Bg decays ||l7| . 
Therefore, although there are weak phases from the c or u quark in the loop, in particular, for B^, the 
effect on the decay phase induced by them is neglegiblly small. 

^Note that the phase convention between B'^ and B'^ is fixed as CP\B'^ —\B'^ > when deriving 



eqs. (yD, duD. 

•^One can check by combining Eqs. ( p^ and (|l^) that all freedoms of phase conventions are canceled 
out completely in ^^j^, including the one between B^ and 
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to x=l^ one will have CP violation of order of 10 ^ for B'^ and 10 ^ for which are 
unobservably small. 

In the models where Eq. (0) is valid, defining ^ = ^ = and using Eqs. 

(I), (p!4|), the time dependent CP asymmetries Eqs. (H) and (H) and time integrated CP 
asymmetries Eq. can be greatly simplified 

41 sm{Amt)Im{0 

^^^^ cos2(iAmt) + |e|2szn2(lAmt) ^ ^ 

.2 gm(Amt)/m(0 

^^^^ |e|2cos2(iAmt) + szn2(lAmt) ^ ^ 



(2 + X2)+X2|e|2 



(17) 



A^. ^ -^-^isfc. (1^ 



(2 + x2)iep + x2 



where = —^{q = d,s for i?^ and respectively). As expected, they are nonzero in 
the presence of CP violating phases. 

We have discussed the CP asymmetries assuming that B^ or B^ mesons are tagged 
before the decay 5^ l^l^{q = d, s) happen f\ Likewise one may also define CP asym- 
metries of the opposite tagging order |[n| which turn out to be just of the opposite sign 
of those defined above, (|1^ and ([T8|). (Eqs. ([15|) and (0) hold for either tagging order.) 
The CP asymmetries not requiring measurement of the time order as one may naively 
imagine to define, however, turn out to be zero because of the the relation Eq. (^ and 
the approximation Eq. ([l^) we have used in our discussions. 

From Eq. (p!7| ) and (0) one can simply get the maximal limit of the CP violating 
observables 

Alil{B'^)U^ = ^=L= (19) 
2 + X2 

which is about 63% for q=d and 5% for q=s. For B^ we know that Xg is experimentally 
larger than 15.7(90% CL)0, so we can neglect the number 2 in the formula and get 

AcpiBg) = 'x^sl' + X^ = 

The situation is clearly quite different for B^ because X^ is just about 0.73. The two CP 
asymmetries for B^ do not exhibit strong correlation. 

In Fig. 1 the correlation between the CP asymmetries of B^ and is plotted scanning 
the parameter space of Cq^ and Cq^ with |Cq. | > 0.1 {i = 1, 2). The points in the figure 
are plotted satisfying the constraints Eq. (|1|). One sees that they do not exhibit strong 

According to the box diagram calculation in SM, the deviation of x from 1 is ^ 10^''(10^'^) for 
Bd(-Bs)[0- So 10~^ is a conservative estimate. 

^An analysis of tagging has been carried out in Rcfs. JTs], |ll]. 
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correlation in the parameter space which is actually implied by the fact that in the most 
of the parameter space, (of order one) is not important at all because of the very 
large X^, while \^d\^ would compete with Xj in the formula Eq. ([T7|). 

We now discuss CP violation of the process in a general 2HDM and MSSM. It has been 
shown [O] that the contributions to the mixing of and from 2HDM or MSSM can be 
significant when the charged Higgs boson mass and tan/3 are small (m/^± < 200 GeV and 
tan/3 ~ 2) or the gluino mass and the squark mass are small (around 100 GeV and 200 GeV 
respectively) and tan/5 is also small. While all other contributions suppressed in the large 
tan/3 limit, the only contribution surviving in this limit is the contribution coming from 
exchanging neutrilino and down-type squarks and the contribution can become important 
only in a very narrow region of down-type squark mass in the low mass spectrum case [19]. 
In the following we limit ourself to discuss CP violation for and B^ decays far away 
from these regions, i.e., in the regions with large tan/? and relatively heavier down- type 
squark mass. Therefore, to a good approximation we can take the mixing to be that 
in SM, i.e., Eq.(|l^. Thus we can use Eqs. ([T5|), (0), (p^) and ( p!8D in the numerical 



analysis. The explicit expressions of the Wilson coefficients Cio, Cq^,Cq2 in a CP softly 



broken 2HDM and MSSM can be found in Refs.pO 



For a CP softly broken 2IIDM[^, the CP violation is depicted by the phase of vacuum 
(i-e., ^ in Ref. [|20[). In Fig. 2 we give the plots of A^p for B^ as functions of varying 
between [0,27r]. Other parameters describing the model are chosen as Mhi = 120 GeV, 
Mh2 = Mh± = 200 GeV, tan/5 = 50 for which the experimental constraints of K — K 
and B — B mixing , r(6 57), r{b ctVt) and are well satisfied. The constraints of 
electric dipole moments (EDMs) of the electron and neutron are also satisfied except for 
in = f P3- ^^^^ ^-P asymmetry can be as large as 20% in vast of the 



range of C,h and can even reach 50%. For B^, the dependence of the CP asymmetry on 
C,H is similar to that for B^ and the CP asymmetry can reach 3%. 

For generic SUSY models, the constraints from the EDMs of the electron and neutron 



on the CP violating phases have been analyzed by many authors|21, 10|. The scenario 



with large tan/?, which we are interested in here, have been discussed in our previous 
papers[|l^]. The constraint of B ^ X^j has also been presented there. In the case of 
low mass spectrum (the lighter stop of order 200 GeV and chargino masses less than 200 
GeV), Cqi and Cq2 are constrained by the B Xg'j decay, because Cq. and C7 (the 
branching ratio of 6 ^ 57 is determined by ) both receive most important SUSY 

contributions from exchanging top squark. An interesting case happens when the SM 
contribution to Cj is completely canceled by the real part of SUSY contributions and 
a considerable imaginary part is left0] ( so that the constraint on C7 is satisfied ) if 
tan/3 is large enough (say, > 30). Cqi and Cq2, in this case, exhibit phases about ±7r/4 
and consequently the absolute value of CP asymmetries for B^ can be significantly larger 
than 30%. CP asymmetries for B^ can also be ±3% in this case. As pointed out in 



Ref. ||10||, the above areas of parameter space are allowed by the EDM constraints due to 
the cancellation among the various contributions to EDMs. For the case of high mass 
spectrum where the B ^ Xs'~f constraint can be safely satisfied and the CP violating 
phases of trilinear term. At, and fi can survive in almost all of their parameter space 
after satisfying the constraints of electron and neutron EDMs, the magnitudes of Cqi 
and Cq2 are also suppressed by the mass spectrum and CP asymmetries can exhibit the 
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correlation depicted in Fig. 1. But for this scenario the branching ratio of the decay would 
not be enhanced large enough, so it is less interesting. In the supergravity(SUGRA) model 
there is another feature which would have important phenomenological implications, i.e., 
because electroweak (EW) symmetry is broken spontaneously by the radiative breaking 
mechanism the masses of the two heavier neutral Higgs bosons are of the same order. 
Hence in general there is a large cancellation happened in the numerator of Eq. ([T3|) in 
SUGRA models. The consequence of it is that for A^^p is greatly suppressed (see Eq. 
(p^)) even in the case of low mass spectrum and the two CP asymmetries for are both 
small (< 10-2). 

With the branching ratios Br{B^ ^l^l) ^'^'^ Br[B^ —>■ given respectively by 



Cbo X 

1 

and 
where 



4m^)\CQ,\^ + \Cq, + 2miCg\^ - 2^1 - Amf[C*Q, x (Cq22 + mC^) + h.c. 



(20) 



:i - Amf)\CQ,\^ + \Cq, + 2miC,\^ + 2^1 - Amf[C*Q, x {Cq^ + 2miC,) + h.c. 



647r3 



Cb; = ^^2F<»-b./I.|a.iVi-4a?. 



we calculate the events A*"* (i=l,2) needed for observing A^p at la in the areas of pa- 
rameter space in which and the branching ratios both have large values and all 
experimental constraints are satisfied. For l=/i, they are order of 10^ and 10^ for B'^ 
and Bg respectively in 2HDM with CP violating phase and tan/?=50 or in SUSY with 
tan/3=30 as well as sparticle masses in a reasonable range. Therefore, 10^° (10^^) B^ (Bg) 
per year, which is in the designed range in the future B factors with 10^- 10^^ B hadrons 
per year [Q, is needed in order to observe the CP asymmetry in i? ^ /i^/i" with good 
accuracy. For l=r, the events A"* are order of 10^ and 10^ for B^ and respectively 
in 2HDM with CP violating phase and tan/5=50 or in SUSY with tan/?=30 as well as 
sparticle masses in a reasonable range. Assuming a total of 5 x 10*^(10^) BdBd (BsBg) 
decays, one can expect to observe ~ 100 identified Bg — > t~^t~ events, permitting a test 
of the predicted CP asymmetry with good accuracy. 

As discussed above, we need to seperate the final state l^lj^ from Ip^lp in order to 
measure CP asymmetry. For l=r, the polarization analysis is straightforward. However, 
detecting tau's is difficult experimentally. For l=/i, in principle one can separate the 
final state UlI^'l from fJ^^fiji by measuring the energy spectra of the electron from muon 
decay[^]. A hl will decay to an energetic ei, which must go forward to carry the muon 



spin, and a less energetic pair of neutrino and antineutrino because the electron is always 
left-handedQ and the energy-momentum and angular momentum are conserved. Due to 
the same reason, for fifi, the relative energies of electron and a pair of neutrino and 
antineutrino are roughly reversed. Therefore, the energy spectra of the electron from 

''In the present case it is quite a good approximation to ignore the mass of electron. 
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muon decay is a powerful /i spin analyzer. However, in practice muons never decay in 
a 4n detector because the lifetime of a muon is long ( cr=659 m). A possible way to 
measure a polarized muon decay is to build special detectors which can make muons lose 
its energy but keep polarization so that the polarized muon decays can be measured. 

In summary, we have analyzed the CP violation in decays — > /"*"/" (q=d,s). While 
there is no direct CP violation, there might be mixed CP violation in the process 

B^^B^^f vs. B^^B^^f. (21) 

It is shown that in the approximation of 1 1 1 =1 the CP violation vanishes in SM. If including 
the correction of order of 10~^ to 1^1 = 1, CP violation is smaller than or equal to 0(10^'^) 
which is unobservably small. In a 2HDM with CP violating phases and MSSM the CP 
asymmetries depend on the parameters of models and can be as large as 40% for B^ and 
3% for Bg. Therefore, an observation of CP asymmetry in the decays B^ — >■ l^l~{q = d, s) 
would unambiguously signal the existence of new physics. 
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Figure Captions 

Fig. 1 The correlation between CP asymmetries for B'^ and B^. 
Fig. 2 A(.p for versus the CP violating phase in 2HDM. 



8 



0.1 




.Q -| I I I I I I I I I 

-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

A'cp(B°d) 



0.8 



0.6 



0.4 



0.2 



i-i 



T3 

O 

^ 
o 

CO 

< 



-0.2 



-0.4 



-0.6 



-0.8 



